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We report on an exhaustive investigation of the dynamical dimer-dimer correlations in imaginary 
time for the quantum dimer model on the triangular lattice using the Green's function Monte 
Carlo method. We show in particular that soft modes develop upon reducing the dimer-dimer 
repulsion, indicating the presence of a second-order phase transition into an ordered phase with 
broken translational symmetry. We further investigate the nature of this ordered phase, for which 
a 12-site unit cell has been previously proposed, with the surprising result that significant Bragg 
peaks are only present at two of the three high-symmetry points consistent with this unit cell. We 
attribute the absence of a detectable peak to its small magnitude due to the nearly uniform internal 
structure of the 12-site crystal cell. 
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I. INTRODUCTION 

Resonating valence-bond (RVB) states constitute a 
major theme in strongly correlated systems, both in the 
context of Mott insulators and of superconductors i£ 
In his milestone paper, Anderson proposed that high- 
temperature superconductors are intimately related to 
Mott insulators and that the pairing mechanism is due 
to "preformed" pairs already present in the strongly cor- 
related insulating phase. 1 Originally, the RVB phase was 
defined by a fully-projected BCS wave function where 
the electron pairs have an arbitrary range. A simplifi- 
cation of the RVB construction was further proposed by 
Rokhsar and Kivclson, who considered an effective quan- 
tum dimer model (QDM) with only local processes and 
orthogonal dimer coverings^ The conditions under which 
such a Hamiltonian could be an accurate description of 
an SU(2) Heisenberg model are not fully understood yet, 
but it is expected to be a reasonable approximation if 
the dimer coverings constitute a natural variational basis, 
and specific quantum dimer models have recently been 
derived for the trimerizcd Kagome antiferromagneli^ by 
Zhitomirsky, 5 and for a spin-orbital model by Vernay and 
collaborators^ 

Regardless of its actual validity for microscopic spin 
models, the QDM of Rokhsar and Kivelson has attracted 
a lot of attention as a promising way to investigate RVB 
physics. In particular, the QDM on the triangular lat- 
tice has been shown to possess a liquid phase with an 
exponential decay of all correlations^ The Hamiltonian 
of that model is 

H = -i£(l//X~ _ l + ^-) 

where the sum runs over all plaqucttes including the three 
possible orientations. The kinetic term controlled by the 



amplitude t changes the dimer covering of every flippable 
plaquette, i.e., of every plaqucttc containing two dimcrs 
facing each other, while the potential term controlled by 
the interaction V describes a repulsion (V > 0) or an 
attraction (V < 0) between dimers facing each other. 

Using the Green's function Monte Carlo (GFMC) algo- 
rithm to probe the ground-state properties of large clus- 
ters, we recently proved the existence, in the thermody- 
namic limit, of topological degeneracy in the disordered 
(i.e., RVB) phased This degeneracy persists over an 
extended parameter range below the Rokhsar-Kivelson 
point V/t = 1. With decreasing V/t, the RVB phase is 
replaced by a crystalline phase with a 12-site unit cell, 
the so-called vT2 x \/l2 phase. Upon further decreas- 
ing the dimer-dimer repulsion, this phase is ultimately 
followed by a columnar phase for sufficiently large and 
negative V/t. 

From the point of view of RVB physics, one of the most 
interesting questions is the nature of the crystallization 
transition between the RVB and the \/l2 x VT2 phases. 
In the present work, we analyze this phase transition by 
extending the GFMC algorithm to calculate dynamical 
dimer-dimer correlations. This has allowed us to extract 
the energy of the first excited state throughout the Bril- 
louin zone (BZ). Having in mind that the \/l2 x \/l2 
phase has dimer order in the thermodynamic limit, the 
dimer gap must close at some k-points of the BZ upon 
entering that phase. An analysis of the gap behavior 
(on the disordered side) and of the structure factor (on 
the ordered side) near the transition point gives a strong 
evidence in favor of the second order phase transition. 

The paper is organized as follows: In Section ^ we 
give a brief explanation of how we calculate the dynami- 
cal correlations and discuss the importance of the choice 
of the finite-size clusters. Then, in Section ITTT1 we discuss 
the excitation spectrum, as well as the location and the 
nature of the phase transition between the RVB liquid 
and the crystalline phase. Finally, Section IIVI is dedi- 
cated to the internal structure of the Vl2 x \/T2 phase. 
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Concluding remarks are given in Section |V| 



II. THE NUMERICAL METHOD 

The numerical investigations are based on a zero- 
temperature Monte Carlo method that filters out the 
high-energy components of a starting wave function | ty G ) 
by the iterative application of G = (A — H), where A is 
a diagonal matrix proportional to the identity. This al- 
lows us to sample directly the ground state l^o)- Indeed, 
in the case of the QDM defined by Eq. all the off- 
diagonal matrix elements are non-positive and, therefore, 
there is no sign problem and the GFMC is numerically 
exacts Moreover, besides the ground-state energy and 
static correlation functions, it is also possible within the 
same method to compute dynamical correlations in imag- 
inary time, such as the dimer-dimer correlation function 
defined by 



D(k,r) 



(y G \d(-k)e- HT d(k)\y ) 



m G \e 



-Ht 



(2) 



where d{k) is the Fourier transform of the dimer operator 
d(r) in real space that is equal to 1 if there is a dimer 
at site r in the (1,0) direction and otherwise. The ex- 
pression of Eq. J2J) can be easily computed by GFMC 
by applying the so-called forward-walking technique^ 
This method to calculate dynamic correlation functions 
has been recently used to investigate an XY spin model 
with ring-exchange^i In practice, after an equilibration 
time, the dimer configurations \x) are statistically sam- 
pled. Then the (diagonal) dimer operator d(k) is applied 
to | re) and the system is evolved for L further steps by the 
forward walking. Within this technique, it is very impor- 
tant to analytically perform the limit of A — > oo^ which 
allows us to speed up the algorithm and define the imag- 
inary time r. In more detail, in the limit A — ► oo, the 
probability of the diagonal moves grows and, therefore, it 
is better to sample directly the number of diagonal moves 
before changing the configuration \x). Note that, in or- 
der to have a finite probability for non-diagonal moves 
for A — > oo, we have to apply G to a power proportional 
to A, i.e., G 0A . Then 



lim (A-H) 



f3A 



-H0 



A— >oo 



(3) 



This algorithm allows us to perform at once all the diag- 
onal moves and then perform directly the non-diagonal 
one, thus saving a huge computational effort. Then 



D{k,Lf3) 



lim 

M—KXS 



{^ G \d{-k)e- HLf3 d{k)e- HMfi \^ G ) 
(^ G \e- H ( L + M ^\^ G ) 



(4) 

and the imaginary time of Eq. (j2J) is given by r = L[3. 
It is easy to demonstrate that for large r, when all the 
high-energy states have been projected out, we have 




D(k, t) oc e 



-A(k)r 



(5) 



FIG. 1: (Color online) Left: Example of cluster used in the 
calculations (here the 12-site one). Right: Brillouin zone of 
the triangular lattice. High-symmetry points are labelled by 
K = (4tt/3, 0), M = (tt, ty/VS), and X = (2tt/3, 0). 



Then, by fitting the large-r behavior of D(k,r), it is 
possible to estimate the gap A(fc) of the first excited state 
at the k-point in the BZ. Note that, within this approach, 
we only find the excitations that have a finite overlap with 
the state d(k)\^ ). In this work, we are only interested in 
the lowest excitation energy, which can be extracted by 
considering only one exponential and fitting the largest-r 
regime. In order to verify the accuracy on A(fc), we have 
also checked the results by fitting with two and three 
exponentials. For extracting more information on the 
spectrum, one may implement more involved methods 
such as the maximum-entropy technique^ which goes 
beyond the scope of the present paper. 

Since we are interested in the transition between the 
disordered RVB phase and the crystalline \/l2 x \/l2 
one, we have decided to work only with clusters defined 
by factorization of the infinite triangular lattice by the 
translations^ 

Ti = iui +lu 2 
T 2 = -hii + 2Zu 2 , 

where u x = (1,0) and u 2 = (1/2, \/3/2) are the uni- 
tary vectors defining the triangular lattice and I is an 
integer. Such clusters can accommodate the proposed 
crystal structure without defects. The number of sites in 
the cluster is N = and to admit dimer coverings I 
must be even. An example of such a cluster for I = 2 is 
shown in Fig.^ The right panel of Fig. Q] shows the BZ 
of the triangular lattice with the high-symmetry points 
K = (4-7r/3,0), M = (7r,7T/V3), and X = (2tt/3,0). 

We remind that the QDM on a toroidal cluster has 
four topological sectors denoted by (0,0), (1,0), (0,1), 
and (1,1): The Hilbert space splits into four subspaces 
that cannot be connected by local dimer flips like the 
ones contained in the Hamiltonian. Given the rotational 
symmetry of the cluster, three of the four topological 
sectors are always exactly degenerate, and the only inde- 
pendent ones are the (0, 0) and the (1, 1), in the notation 
of Ref. 0. The ground state belongs to the former or the 
latter depending on the parity of 1/2 A It is worthwhile to 
mention that, within GFMC, we can fix the topological 
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FIG. 2: (Color online) Typical results for the dynamical 
dimer-dimer correlations D(k, r) for the 432-site cluster at 
X = (2-7T-/3, 0) (upper panel) and K = (47r/3, 0) (lower panel) 
for a few selected values of V/t. 



sector during the simulation and, therefore, we can study 
the different sectors independently. 

In the following, we will present the results for the 
dimer gap along particular symmetry directions of the 
BZ, and in particular the size scaling of the k-points 
where soft modes are expected to develop upon enter- 
ing the crystalline Vl2 x \/l2 phase. 



III. THE RVB - CRYSTAL TRANSITION 

In their original work, Moessner and Sondhii estimated 
that the transition from the disordered to a crystalline 
phase occurred around V/t ~ 0.7. A later analysis of 
the topological splitting on finite-size clusters suggested 
a slightly larger value (0.7 < V/t < 0.85)£ The analy- 
sis of dynamic correlations reported in our present work 
provides a more accurate estimate, V/t = 0.82 ± 0.03, as 
will be demonstrated below. 

In the crystalline phase, the existence of a broken 
translational symmetry naturally implies a folding of the 
BZ and soft modes at particular k-points. In particu- 
lar, for the 12-site unit cell proposed by Moessner and 
Sondhi, we expect to find soft modes at the X, M, and 
K points (see Fig.^l. So we now turn to an investigation 
of the excitation spectrum as revealed by our GFMC ap- 
proach, with emphasis on the behavior around these high 
symmetry points. 

We perform a systematic analysis of the dynamical 
correlations for clusters with up to 432 sites and for 
0.55 < V/t < 1.0. We report in Fig. [3 the results for 
a 432-site cluster (JV = 3l 2 with I = 12) at X and K 
for different values of V/t. Although in principle the 
long-time dynamics should always give access to the en- 
ergy of the first excited states, in practice the error bars 
grow with time, and two types of data sets may be dis- 
tinguished. In the first, favorable case, an exponential 
decay governed by a single exponent shows up long be- 
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FIG. 3: (Color online) Lowest-energy excitation spectra for 
the 432-site cluster for different values of V/t. 



fore the error bars become significant, and the gap can be 
extracted very accurately (see upper panel of Fig. 0) • In 
a logarithmic scale, this shows up as a linear behavior of 
D(k, r) prior to the development of large error bars. This 
is the case throughout the BZ at the Rokhsar-Kivelson 
point (V/t = 1), and everywhere except around the K 
point when V/t < 1. In the other case, the error bars 
start to grow before a simple exponential decay has set 
in. In logarithmic scale, this shows up as a persistent 
curvature. An example is reported in the lower panel of 
Fig. |21 The absence of a clear exponential regime pre- 
vents us from extracting a reliable value for the energy 
of the first excited state. 

The results of this analysis are depicted in Fig. [3 where 
we show the excitation spectrum of the 432-site cluster- 
in the topological sector (0, 0) (the one that contains the 
ground state) for different values of V/t. The missing 
points around K correspond to cases where the extrac- 
tion of the energy was impossible. This spectrum calls 
for several comments. 

First of all, the spectrum is fully gapped close to the 
Rokhsar-Kivelson point, in agreement with the results 
obtained previously by Ivanov at the Rokhsar-Kivelson 
point^ and consistent with a liquid ground state with 
no long-range order. 

Secondly, upon decreasing V/t, two soft modes clearly 
appear at X and M. These two k-points are compati- 
ble with the 12-site unit cell structure of the VT2 x \/l2 
phase, proposed by Moessner and Sondhi^ It turns out 
that the numerical results concerning X and M are very 
accurate and allow us to obtain a rather precise deter- 
mination of the transition point. In Fig. 01 we show the 
size scaling of these two gaps. At the Rokhsar-Kivelson 
point, the curves saturate quickly to a finite value in the 
thermodynamic limit. This value becomes smaller by de- 
creasing the ratio V/t, and we can estimate that the tran- 
sition occurs around V/t ~ 0.8, not far from our previous 
estimate based on the topological gap£ This constitutes 
one of the important results of this paper: The transition 
out of the RVB phase is second order, as revealed by the 
development of soft modes upon entering the crystalline 
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FIG. 4: (Color online) Size scalings of the gap at the X 
(upper panel) and M (lower panel) points for different values 
of V/t. 
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FIG. 5: (Color online) Static dimer-dimer correlation at 
V/t = 0.55 for different cluster sizes (upper panel). Size scal- 
ing of the dimer order parameter square D\ (calculated at the 
X point) for different values of V/t (lower panel). 



phase. 

Thirdly, and quite surprisingly, we could not get any 
evidence for the emergence of a soft mode at the K point, 
the third high-symmetry point compatible with the 12- 
site unit cell. On the contrary, the slope of D(K,t) (in 
logarithmic scale) at the largest accessible times is still 
very large, so that there is no evidence whatsoever of 
a small energy scale at that point. We can think of two 
possibilities to explain this result: i) The overlap between 
the first excitation and the starting state d(fc)|^'o) is so 
small that this excitation does not give any detectable 
contribution to the dynamics in the accessible time range; 
ii) The symmetry of the ordered phase is such that no soft 
mode develops at the K point. To clarify this point, we 
now turn to a detailed investigation of this crystalline 
phase. 



IV. THE 712 x 712 PHASE 

The unexpected behavior found at K is actually not in 
contradiction with any of the previous numerical studies 
of the static correlations in that phase: Moessner and 
Sondhi detected a peak at X and noticed it was compat- 
ible with a 12-site unit cell, but they did not calculate 
the full k-dependence of the static dimer-dimer correla- 
tionsi In our previous paper^ we concentrated on the 
RVB phase and showed that the correlations were indeed 
decreasing very fast in that phase, but we did not inves- 
tigate in details the correlations in the ordered phase. So 
the first thing to do is to go back to the static correla- 
tions and to calculate the full k-dependence of the static 
dimer-dimer correlation function defined by: 



(V \d{-k)d{k)\* ) 
(*o|*o) 



(6) 



The results for different cluster sizes and V/t = 0.55 
are depicted in Fig. [S] Clearly, these results confirm the 
anomaly of the K point: While strong peaks are present 
at X and M and grow with the size of the cluster, no peak 
could be detected at K. Notice that for this kind of static 
correlations, there is no fitting procedure, in contraxt to 
the extraction of the gap in the dynamical calculations. 

Before trying to understand what is going on at the K 
point, let us look more closely at the X and M points. 
For small enough V/t, the size of these peaks grows lin- 
early with the number of sites of the cluster, a clear in- 
dication that Bragg peaks develop. Quantitatively, we 
extract the square of the order parameter defined by: 



■Pstat(fc) 

N ' 



(7) 



for both X and M from the finite-size scaling of these 
peaks. The results are depicted in Fig. [SJ The extrap- 
olated values are very small even well inside the crys- 
talline phase, but given the accuracy of the data and 
the smoothness of the size scaling, these values are reli- 
able, and indeed they again point to a second-order phase 
transition located between V/t ~ 0.8 and 0.85 where the 
order parameters continuously vanish at both points (see 
Fig.|HJl. So the existence of a crystalline phase with Bragg 
peaks at X and M is well established. 

A similar analysis of the Bragg peak at the K point 
gives a zero value for D\, within the error bars. For 
example, well inside the crystal phase at V/t = 0.55 we 
can estimate an upper bound D\ < 1.5 x 10~ 6 . 

As mentioned at the end of the previous section, there 
are two possible explanations for the absence of the K 
peak. The first possibility is an exact cancellation of the 
K peak due to the symmetry of the crystal phase. The 
second one is just a strong suppression of the K peak so 
that its magnitude is below our numerical precision. 

The symmetry cancellation of the K peak might occur 
in a crystal phase with a three-fold rotational symmetry 
about a lattice siteil However this type of crystal sym- 
metry seems to be unlikely, since it cannot accommodate 
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FIG. 6: (Color online) V/t dependence of the thermodynamic 
values of D\ and of the gap at the X and M points. 



any individual dimer or plaquette covering with the sym- 
metry of the crystal (due to the single site in the center 
of the three-fold symmetry) . Note that the VT2 x \f\2 
crystal structure proposed by Mocssner and Sondhii^ has 
a three-fold symmetry about a triangular plaquette cen- 
ter and should not lead to an exact cancellation of the 
K peak. To get a further insight in the symmetry of the 
\/T2 x VT2 phase, we have calculated the static density 
profile in the presence of a local perturbation breaking all 
symmetries of the lattice. In this way, instead of obtain- 
ing a uniform linear combination of all symmetry related 
crystal states, the system selects one (or sometimes sev- 
eral) crystal state(s) favored by the perturbation, giving 
access to a real-space snapshot of the broken-symmetry 
ground state far enough from the perturbation. A typical 
result obtained by considering local potentials for three 
bonds related by a 27r/3 rotation around a site that favor 
or disfavor the presence of a dimer is depicted in Fig. [7| 
The 12-site diamond of the type proposed by Mocssner 
and Sondhii& is clearly visible. Other choices of potential 
perturbation sometimes produce less symmetric periodic 
patterns which we interpret as superpositions of several 
crystal states. However, in no pattern were we able to 
detect a three-fold symmetry about a lattice site. There- 
fore we conclude that the symmetry-cancellation scenario 
is not realized and that we should interpret the absence 
of a visible K peak as a strong suppression of the cor- 
responding Fourier harmonic in the \JY1 x ^fY2, state of 
Mocssner and Sondhiii& 

Interestingly enough, while the unit cell is clearly visi- 
ble in Fig.[?l there is no internal modulation of the dimer 
density inside the resonating 12-sitc hexagons. Motivated 
by these observations, we have performed a calculation 
with a trial wave function containing a uniform RVB 
superposition of all dimer configurations within 12-sitc 
hexagons and no dimers between the hexagons. The re- 
sulting peak weights D 2 X and D 2 M are of the same order 
of magnitude, while the peak weight D\ is more than 
two orders of magnitude smaller. Such a huge difference 
between the peak weights would explain our numerical 
findings about the absence of the K peak for both dy- 
namic and static correlations. For illustration, the static 




FIG. 7: (Color online) Mean value of the dimer operator for 
the 192-site cluster for V/t — 0.6 in the presence of the lo- 
cal perturbation (shown in the inset) placed at the rightest 
site. Dashed lines correspond to correlations lower than the 
uniform distribution (i.e., 1/6) and continuous lines to corre- 
lations larger than 1/6. The thickness of the bonds indicates 
how different the amDlitude is from the uniform distribution. 
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FIG. 8: (Color online) Static dimer correlations in the trial 
wave function constructed as the product of uniform RVB 
superpositions of dimer covering of 12-site hexagons. 



dimer-dimer correlations for our trial wave function are 
shown in Fig. |SJ They qualitatively resemble the actual 
correlations in the crystal phase of Fig. [5] with no visible 
K peak. 



V. DISCUSSION 

To summarize, we have studied the phase transition 
between the RVB and the crystal states of the QDM 
model on the triangular lattice by analyzing the dimer 
excitation spectrum with the help of the Green's func- 
tion Monte Carlo method. First of all, we have shown 
that upon approaching the transition between the RVB 
phase and the ordered phase, two soft modes develop in 
the ordered phase, while two Bragg peaks continuously 
disappear in the RVB phase. The continuous behavior 
of the dimer gap and of the Bragg peaks at the transi- 
tion point strongly suggests the second order of the phase 
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transition, consistent with the theory of Moessner and 
Sondhi. 16 Furthermore, our analysis provides a more ac- 
curate location of the transition point V/t = 0.82 ± 0.03 
than in previous studies. 

A peculiar feature of the crystal phase is the apparent 
absence of a singularity (Bragg peak or soft mode) at 
the zone corner. We interpret this effect as a strong re- 
duction of this Fourier harmonics in the crystal without 
much short-length structure in the large 12-site crystal 
cell. In other words, the \/l2 x \/l2 crystal may be bet- 
ter approximated by 12-site hexagons without internal 
structure, and not by resonating 4-site plaquettes. 

It is interesting to compare these new findings with 
the phase diagram of the same model on the square lat- 
tice. In that case, the Rokhsar-Kivelson state realized at 
the Rokhsar-Kivelson point (V/t = 1) is immediately re- 
placed by a plaquette phase, until columnar dimer order 
is finally stabilized for small enough repulsion. The pla- 
quette state is a local RVB state on a square plaquette, 
and in that sense, it is similar to the local RVB state on 
the 12-site diamond. But the analogy can be extended 
further. The 4-site plaquette is the smallest subunit that 
allows resonance between pairs of dimers and that keep 
the fundamental point group symmetry of the square lat- 
tice, the C4 axis. The same requirements naturally lead 
to the 12-site diamond cluster for the triangular lattice: 



A 4-site plaquette does not have C3 symmetry, and the 6- 
site plaquette that has the shape of a large triangle does 
not allow pairs of dimers facing each other. So one can 
conjecture that there is a general tendency for such mod- 
els, when the dimer repulsion is no longer strong enough 
to stabilize the fully resonating state, to develop local 
RVB order in symmetric unit cells that allow sufficient 
resonance. 

Finally, although wc know quite a lot by now about 
the QDM on the triangular lattice, a number of very 
interesting questions remain, among which the role of 
temperature, the nature of the excitation spectrum in 
the RVB phase, and the possible role of visons at the 
transition to the ordered phase. 
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